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摘 ， 要 :在 强 外 界 载荷 下 ,在 海洋 平台 - 立 管 多 体系 统 中 的 海洋 平台 和 立 管 间 的 相互 耦合 作用 会 加 
强 。 从 而 导致 整个 系统 的 非 线性 增强 。 考 虑 到 动力 定位 海洋 平台 - 立 管 多 体系 统 的 强 非 线性 ,结合 
真实 的 海上 施工 工程 背景 和 凝集 质量 法 ,基于 Python 镶 谋 编程 和 OrcaFlexAPI 模块 的 组 成 以 及 运用 
规则 ,对 OrcaFlex 进行 了 局 部 的 二 次 开发 ,建立 了 一 种 基于 无 迹 卡尔 受 滤波 的 UKF-PID 控制 的 动力 
定位 平台 - 立 管 多 体 耦 合 系统 。 最 终 建 立 了 无 迹 卡 尔 受 滤波 模式 下 PID 控制 的 动力 定位 海洋 平台 - 
5 立 管 刚 柔 多 体 模 型 ,并 对 该 模型 在 特定 海洋 环境 下 进行 了 动态 仿真 。 计 算 结果 对 于 有 具体 工程 实践 
有 着 一 定 的 指导 意义 。 
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Abstract :Under strong external loads, the intercoupling effect between the offshore platform and the riser 
in the offshore platform-riser multi-body system will be significantly enhanced. This leads to the non-line- 
arity enhancement of the whole system. With the consideration of the strong non-linearity of the dynamic 
positioning (DP) offshore platform-riser multi-body system, combined with the actual offshore construction 
engineering background and the lumped mass method, based on Python Language embedded program- 
ming, with the basis of the operation principle of OrcaFlex API and the composition of its modules, the 
partial secondary development of OrcaFlex is carried out, and a UKF-PID control method based on Un- 
scented Kalman filter is proposed to realize the control of the dynamic positioning offshore platform-riser 


multi-body coupling system. Finally, the model of PID-controlled dynamic positioning offshore platform- 
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riser rigid-flexible multi-body system under unscented Kalman filter mode is established and the dynamic 


simulation is carried out under specific marine environment. The calculations obtained have certain guid- 


ing significance for specific engineering practice. 


Key words :dynamic response; offshore platform-riser; Python ; unscented Kalman filter; OrcaFlex; lumped 


mass method 


现代 海洋 工程 在 具体 工程 实践 中 所 面 对 的 技术 
问题 越 来 越 复杂 "i 。 作 为 一 种 常见 的 采油 系统 , 动 
力 定位 海洋 平台 - 立 管 系统 的 总 体 动力 学 特性 是 一 
个 值得 研究 的 课题 “1 。 

相对 于 当前 的 动力 定位 系统 来 说 ,传统 的 卡尔 
曼 滤波 以 及 扩展 卡尔 曼 滤波 在 海洋 平台 动力 定位 系 
统 中 的 应 用 非常 广泛 "1。 传 统 的 卡尔 曼 滤 波 
(Kalman filter, KF) 扩展 卡尔 曼 滤 波 (extanded Kal- 
ma fliter,EKF ) 都 是 将 问题 转化 为 线性 高 斯 模型 ,从 
而 天 便 对 问题 进行 求解 运算 。 但 对 于 非 线性 问题 ， 
双 E 是 KF 还 是 EKF ,都 存在 计算 量 较 大 、 易 受到 线 
性 涡 差 影响 的 问题 ,而 线性 误差 的 引入 会 降低 模型 
的 芋 度 。 因 此 ,引入 无 迹 卡 尔 曼 滤波 (UKF) , 它 不 同 
了 RDKF ,由 于 没有 将 高 阶 项 忽略 ,故而 这 种 方法 对 统 
做 的 近似 精度 较 高 。UKF 的 效果 理论 上 可 以 达到 
二 阶 EKFCEKTF 将 非 线性 模型 在 状态 估计 值 附近 做 
秦 盎 级 数 展开 ,并 在 一 阶 稚 断 ) 才能 达到 的 效果 *”。 


1 动力 定位 海洋 平台 - 立 管 系统 
Fig.1 Dynamic positioning offshore platform-riser system 

在 外 界 载荷 变化 较为 剧烈 时 ,对 于 海洋 平台 - 立 
管 这 种 刚 柔 多 体系 统 来 说 ,很 容易 出 现 非 线性 结果 
线性 化 的 现象 。 为 克服 这 一 缺点 ,将 基于 经 典 的 
PID 控制 方法 '" ,结合 无 迹 卡 尔 曼 滤波 的 优点 ,以 
Python 作为 镶嵌 语言 ,基于 OrcaFlex 设计 一 种 UKF- 
PID 动力 定位 海洋 平台 - 立 管 多 体 耦 合 系统 ,并 对 其 
进行 总 体 分 析 。 


1 海洋 平台 运动 的 控制 


1.1 海洋 平台 运动 的 空间 模型 与 无 迹 卡 尔 曼 滤波 
对 于 动力 学 问题 而 言 ,常见 的 问题 是 线性 问题 


和 非 线 性 问题 。 线 性 系统 和 非 线 性 系统 最 明显 的 区 
别 方法 就 是 线性 系统 遵从 垣 加 原理 , 而 非 线性 系统 
不 然 。 一 般 来 说 ,振动 的 系统 总 是 非 线 性 的 。 其 中 ， 
非 线 性 又 包括 弱 非 线性 和 强 非 线性 。 强 、 弱 非 线 性 
要 根据 具体 的 研究 对 象 和 研究 目的 来 区 分 。 关 于 强 
弱 非 线性 的 定义 ,可 用 如 下 微分 方程 控制 的 非 线性 
振动 方程 来 进行 区 别 ""。 
dx  ， ) 
a i Wh) (1) 
当 式 (1) 中 的 参数 s 是 较 大 的 数 时 ,该 方程 就 
是 一 个 强 非 线性 方程 。 而 海洋 平台 - 立 管 刚 柔 多 体 
系统 以 及 该 系统 所 处 的 海洋 环境 ,决定 了 该 系统 在 
运行 时 会 出 现 强 非 线 性 现象 。 
描述 动力 定位 海洋 平台 的 相关 运动 需要 定义 两 
个 坐标 系 。 一 个 是 地 球 惯性 坐标 系 , 也 称 总 体 坐标 
系 ,主要 用 来 描述 海洋 平台 的 位 置 和 出 向 向 量 x = 
[x,y,W] ; 男 一 个 是 固定 在 海洋 平台 上 和 海洋 平台 


一 起 运动 的 海洋 平台 局 部 坐标 系 ,其 坐标 原点 常设 
置 在 海洋 平台 重心 ,用 来 描述 海洋 平台 的 速度 向 量 
p=[u,v,r] tl, 
动力 定位 海洋 平台 低频 运动 模型 可 表示 为 
n=R(y)r (2 
My = -Dyv+u+R (yb+w, (3) 
b= -T, b+w, (4) 


式 中 :MM 为 包含 附加 质量 的 海洋 平台 惯性 力矩 阵 ;已 
为 海洋 平台 的 线性 阻尼 和 矩阵 ;uw 为 平台 受到 的 控制 
力矩 阵 ;2 为 描述 低频 环境 干扰 和 未 建 模 动 态 的 偏 
差 项 ;w, 、w, 为 零 均 值 白 噪 声 随机 过 程 。 灵 (水 ) 为 海 
洋 平 台 的 旋转 矩阵 , 且 


cosyy -siny 0 
R(y)=| sinfy cosy | (5) 
0 0 1 


引入 海洋 平台 平行 坐标 系 ,其 坐标 原点 与 地 球 
坐标 系 的 原点 重合 ,坐标 系 各 坐标 轴 与 海洋 平台 局 
部 坐标 系 平行 。 用 力 ,表示 平行 坐标 系 下 海洋 平台 
运动 向 量 , 其 与 地 球 坐标 系 下 运动 量 的 转换 关系 为 
,=R" (小) 和 ,从 而 =v, 进而 得 到 平行 坐标 系 下 海 
洋 平台 的 线性 低频 运动 状态 空间 模型 , 即 
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| =4 Xi +Biu +E, w: (6) 


yr.=H, xL+n 
式 中 : xz = [wy,,v ] 为 状态 变量 ;y1 为 被 控 输出 变 
量 ;wi 为 包含 偏差 项 5 在 内 的 扰动 向 量 ;ni 为 测量 
高 斯 白 噪声 。 式 中 各 和 矩阵 定义 如 下 。 


0, ,3 Ls 0, ,3 
a 

0, -MD M 

03 
EE 0, .3] (7) 


动力 定位 海洋 平台 高 频 运动 状态 的 空间 模型 可 
表示 为 


之 =4 E+E,w; (8) 

n=Cé (9) 

式 电 :6 = [EX ,ys ,Ws] 表示 波浪 运动 的 向 

量 sw; 为 零 均 值 白 噪声 过 程 ;m= [x,,y,,W,] 为 
小江 的 高 频 运动 向 量 。 各 矩阵 定义 如 下 。 

03 .3 Ts 03 

| 

ES C, = [03x1 LT.1] (10) 

式 电 :0 = diag | wo , Wo, woo) ; A = diag { 20w0, 

2 名 名 ,266w01 ,这 里 oo 为 波浪 峰值 频率 ,&; 为 相对 

阻 届 系数 ,一 般 取 值 为 0.05 ~0.2'。 

,综合 考虑 海洋 平台 的 低频 运动 模型 式 (2) ~ 式 

( 怨 以 及 高 频 运动 模型 式 (8) ~ 式 (9) ,得 到 动力 定 

人 i 唾 洋 平台 的 状态 估计 与 滤波 模型 , 即 

I 下 =Ax +Bu +Ew 


SS 4。 一 


(11) 
一 y=Hx+n 
式 惠 :x =[&" ,wy ,b",v"]" 为 系统 的 状态 向 量 ;w = 
[wi ,w，" ,w3 ] 为 系统 的 过 程 噪声 向 量 ;y = + 
ms。+n 为 包含 海洋 平台 的 高 频 和 低频 运动 的 测量 模 
型 ,这 里 n 为 系统 的 测量 噪声 ; 式 中 各 和 矩阵 定义 详 见 
文献 [ 14 ] 。 
基于 动力 定位 船舶 状态 估计 与 滤波 模型 式 
(11) ,将 其 改写 成 时 间 离 散 形 式 , 即 
Xi =Dx, +Au,+T w, 
=f(x,u) +Tw, (12) 
y=Hx, +n,=h(x,) +n, 
式 中 :$=exp(A :At),A=A- (D-DDB,T=A-!: 
(BD-7)E, 这 里 Az 为 采样 时 间 ;xi 为 系统 在 时刻 
的 状态 向 量 ;y; 为 上 时 刻 系统 状态 测量 值 。 系 统 噪 
声 w; 为 零 均 值 高 斯 白 品 声 ,其 方差 为 0 , 量 测 噪 声 
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n; 是 与 系统 噪声 相互 独立 的 高 其 
为 Ri。 

无 迹 卡尔 曼 滤 波 的 核心 在 于 UT 变换 (unscent- 
ed transform ,UT) ,也 称 UT 变换 。 具 体 过 程 为 ,首先 
通过 UT 变换 进行 非 线 性 模型 的 状态 与 误差 协 方差 
的 递 推 与 更 新 ,然后 再 采用 标准 Kalman 滤波 的 方法 
进行 滤波 设计 。 对 于 式 (11) 所 表达 的 动力 定位 非 
线性 离散 系统 ,应 用 无 迹 卡 尔 曼 滤波 进行 状态 估计 
和 滤波 的 过 程 如 下 。 

1) 定 义 状态 和 协 方差 初 值 。 

to=E(xo),Po=E[(xo -Xo) (xo -2) ] (13) 

2) 计算 左 时 刻 sigma 点 并 选取 合适 的 权 值 。 

PE 


白 噪声 ,其 方差 


Xt1— Vn +AVP,i] (14) 
Wo =A/(n+A) 
Ws=A/(n+A)+(1-o +B) (15) 


W”=W=A/[2(n+A)],i=1,2,.…,2n 

式 中 ;WWW(i=0,1,…,2n) 分 别 表 示 sigma 点 均值 
和 方差 的 权 值 ;n 为 系统 的 维 数 ;A = [a (n+xk) - 
nj 为 比例 因子 ,a 决定 了 sigma 点 在 均值 点 周围 的 分 
布 状态 , 常 在 10“ ~1 之 间 取 值 ,kx 通常 取 为 0;B 决 
定 了 先 验 状 态 估 计 值 的 分 布 状 态 , 对 于 高 斯 分 布 , 取 
B=2 为 最 优 。 

3 ) 时 间 更 新 :通过 状态 和 测量 方程 (12 ) 对 sig- 
ma 点 进行 无 迹 变换 ,得 到 先 验 的 状态 和 测量 输出 的 
预测 值 ,及 其 协 方差 的 预测 值 。 

4) 测 量 更 新 :计算 UKF 增益 矩阵 Ki (K, =P … 
(PY ) 一 ) .状态 估计 值 x 及 状态 误差 协 方差 矩阵 Pi 。 

通过 以 上 几 个 步骤 ,在 每 个 采样 时 刻 不 断 重复 
进行 , 且 将 每 个 时 刻 的 UKF 状态 估计 滤波 值 作为 预 
测 未 来 海洋 平台 动态 的 初始 条 件 。 
1.2 ”Python 锐 髓 语言 下 OrcaFlexAPI 中 实现 对 海 
洋 平台 的 UKF-PID 控制 

当 Python 与 OrcaFlex 在 OrcaFlexAPI 中 进行 结 
合 时 ,对 于 每 一 个 新 的 时 间 步 ,要 在 开始 阶段 就 将 
info. NewTimeStep 的 参数 属性 设置 为 True。 输 入 为 
上 一 时 刻 包含 环境 干扰 、 建 模 误 差 造成 的 海洋 平台 
状态 值 ;输出 为 滤波 后 得 到 的 海洋 平台 状态 值 。 再 
结合 无 迹 卡 尔 曼 滤波 ,在 经 典 PID 控制 系统 通过 外 
部 函数 将 DP Vessel 的 横 荡 、 纵 荡 及 船舶 出 摇 与 目标 
值 进行 比较 ,经 过 控制 方程 的 计算 ,得 到 DP Vessel 
所 需 的 回复 反 力 及 反 力 矩 , 根 据 推 力 分 配 相 关 原 则 


进行 推力 分 配 。 


2 动力 定位 海洋 平台 - 立 管 刚 柔 多 体 
系统 的 建立 


2.1 立 管 的 离散 模型 

关于 系统 中 立 管 的 离散 模型 与 动态 响应 的 计算 
方法 ,运用 的 是 凝集 质量 法 。OreaFlex 也 正 是 基于 
凝集 质量 法 对 管线 建立 模型 的 。 关 于 凝集 质量 法 的 
详细 情况 ,可 参考 作者 已 经 发 表 过 的 文献 [15-18] 。 
2.2 动力 定位 海洋 平台 - 立 管 多 体 模型 的 建立 

整 根 立 管 的 总 长 度 为 531 m, 在 对 立 管 进行 单 
元 划分 时 ,如 果 立 管 的 分 段 单元 长 度 过 长 ,就 会 使 得 
立 管 的 刚度 过 大 ;如 果 立 管 的 分 段 单元 长 度 过 短 , 那 
么 就 会 造成 凝集 质量 节点 过 多 , 立 管 就 会 失去 其 自 
身 悉 有 的 刚性 ,造成 力学 特性 失真 。 
< 二 除去 最 上 端的 11 m, 剩 余 的 520 m 为 立 管 主体 。 
为 便于 下 文 的 描述 ,最 上 端的 11 m 称 为 立 管 的 滑 移 
接 买 端 , 滑 移 接头 与 下 面 的 立 管 主体 的 连接 处 称 为 
直 管 的 顶端 , 立 管 与 海 床 的 匀 接 处 称 为 立 管 的 底 端 。 
诡 必 的 分 段 情 况 如 表 1 所 示 。 立 管 的 相关 参数 如 表 
2 隋 示 。 


表 1 立 管 的 单元 划分 


Tab.1 Unit division of riser 


岁 称 长 度 区 间 /m 单元 长 度 /m 单元 数量 
滑 移 接头 0~11 11 1 
内 芒 敏 感 区 11 ~51 2 20 
其 你 部 分 51 ~531 5 96 
表 2 立 管 的 相关 参数 
Tab.2 Parameters of riser 

名 称 大 小 
外 径 /m 0.65 
内 径 /m 0.6 
弯曲 刚度 /(kN .ma ) 滑 移 接头 508 944 
轴 向 刚度 A(kN . m7) 滑 移 接头 1 
扭转 刚度 /(kN : nm ) 滑 移 接头 391 500 
泊 松 比 0.3 
密度 滑 移 接头 0.385 tm 

立 管 主体 7. 85 tm 
杨 氏 模 量 /MPa 立 管 主体 21 200 
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立 管 顶 端 张 力 必须 能 够 保证 立 管 主体 不 会 发 生 
届 曲 。 因 此 ,在 建 模 过 程 中 ,4 个 Link 单元 被 环 向 均 
布 于 立 管 滑 移 接头 端 ,来 对 立 管 顶 端 张力 进行 补偿 
和 缓冲 ,因此 4 个 Link 单元 等 效 于 4 个 张 紧 器 , 且 每 
个 张 紧 器 的 冲程 均 为 +4 m。Link 单元 弹性 刚度 与 
其 被 拉 伸 后 的 长 度 线 性 相关 ,如 表 3 所 示 。 

表 3 Link 单元 的 参数 
Tab.3 Parameters of Link 


阻尼 /A(kKN. m-1.s) ” 拉 伸 后 的 长 度 /m ”弹性 刚度 /A(kN . m-!) 
5 0~4,12~16 50 


3 4~12 25 


海水 密度 为 1025 kg/m ,水 深 为 500 m, 海 床 类 
型 为 水 平 海 床 ;海流 流速 为 0.5 m/s, 流 向 为 270°, 流 
向 与 浪 向 相同 ,选择 JONSWAP 波谱 为 波浪 载荷 计 
算 波 谱 , 波 浪 谱 的 具体 参数 如 表 4 所 示 。 
表 4 JONSWAP 波谱 的 相关 参数 
Tab.4 Parameters of JONSWAP spectrum 


谱 峰 因子 y 


4.7934 


相关 系数 a Ol Or2 H/s 人 了 As 
0.0086 0.07 0.09 水 9 


该 平台 有 8 个 浮 简 ,排水 量 为 10 000 +, 海 洋 平台 
定位 坐标 为 整体 坐标 系 (Target X= +30 m,Target 了 
= -20 m) , 转 出 角 为 90°, 具 体 参数 如 表 5 所 示 。 

表 5 海洋 平台 的 相关 参数 


Tab.5 Parameters of offshore platform 


吃水 /m 转动 惯量 /(t + m?) 
24.38 7, =500 000 人 =7 000 000 1.=7 000 000 


将 海洋 平台 的 位 移 RAO 及 二 阶 波浪 传递 函数 


QTFs 经 过 ANSYS/AQWA 计算 后 导入 海洋 平台 - 立 
管 多 体系 统 模型 中 来 计 及 波浪 载荷 及 二 阶 波浪 漂移 
力 对 系统 运动 的 影响 ” 。 模 型 建成 后 如 图 2 所 示 。 


到 2 UKF-PID 动力 定位 海洋 平台 - 立 管 系统 模型 示意 图 
Fig.2 Model of UKF-PID dynamic positioning 


off shore platform-riser system 
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3 ”模拟 结果 


3.1 均匀 管内 流 时 的 计算 结果 
3.1.1 不 同 管 内 流 流速 作用 下 立 管 的 张力 

图 3 为 不 同 管内 流 作 用 下 的 立 管 张力 响应 图 。 
观察 图 3 发 现 , 当 管内 流 为 均匀 流 时 ,在 管内 流 流 速 
不 大 时 , 随 着 管内 流 流速 的 增加 ,与 单一 PID 控制 下 


所 得 的 计算 结果 相 比 ,加 入 无 迹 卡 尔 曼 滤波 后 并 未 
改变 立 管 张力 大 小 、 在 时 域 上 波动 情况 、 立 管 张力 大 
小 沿 立 管 长 度 分布 的 情况 和 立 管 张力 变化 的 协调 性 
和 同步 性 。 且 可 以 发 现 , 在 两 种 控制 模式 下 , 当 管 内 
流 流速 较 低 时 , 随 着 管内 流 流速 的 增加 , 立 管 张力 的 


3CB2 不 同 管内 流 流速 作用 下 立 管 的 这 短 和 则 率 
加 图 4 为 不 同 管内 流 作 用 下 立 管 的 弯 矩 和 曲率 分 


了 入 无 迹 卡 尔 曼 滤波 后 ,使 得 立 管 整体 沿 长 度 方 
向 淖 受 的 弯 矩 大 大 增加 ,但 并 未 改变 立 管 上 发 生 较 
大 要 曲 的 位 置 。 

(0 当 管内 流 为 均匀 流 时 ,在 管内 流 流速 较 低 时 , 随 
着 管内 流 流速 的 增加 ,与 单一 PID 控制 下 的 计算 结 
果 独 比 , UKF-PID 动力 定位 系统 作用 下 , 立 管 同一 位 


置 琢 弯 和 矩 变化 不 大 。 这 说 明 在 UKF-PID 动力 定位 
系统 作用 下 ,在 管内 流 流 速 较 低 时 , 随 着 管内 流 流 速 
的 增 大 , 立 管 承 受 的 弯 和 矩 变化 不 大 。 因 此 ,UKF-PID 
动力 定位 系统 的 存在 并 没有 改变 立 管 轴 向 张力 的 分 
布 和 变化 情况 ,而 是 改变 了 立 管 承受 横向 弯曲 载荷 
的 分 布 及 变化 情况 。 

当 管内 流 为 均匀 流 时 ,对 比 观察 单一 PID 控制 
与 UKF-PID 动力 定位 系统 作用 下 不 同 管内 流 流速 
下 的 立 管 弯 和 矩 标 准 差 和 曲率 标准 差 沿 长 度 方向 分 布 
情况 发 现 ,无 迹 卡 尔 曼 滤 波 的 加 入 增强 了 立 管 的 曲 
率 变 化 的 非 线 性 和 曲率 变化 在 整个 长 度 方 向 的 非 协 
调 性 。 由 于 立 管 曲率 标准 差 沿 立 管 长 度 阶梯 型 分 布 
的 增强 ,对 于 某 一 管内 流 流速 下 的 立 管 ,其 曲率 变化 
的 同步 性 和 协调 性 沿 着 长 度 方向 分 布 的 分 段 性 也 更 
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加 明显 。 也 就 是 说 ,对 于 曲率 变化 而 言 , 沿 立 管 整 体 
长 度 方向 的 曲率 变化 的 协调 性 与 同步 性 有 所 降低 ， 
但 其 在 某 些 特定 长 度 范围 内 的 曲率 变化 仍 具 有 一 定 
的 协调 性 和 同步 性 ;不 同 的 是 ,曲率 变化 同步 性 较 好 
的 大 长 度 区 间 和 大 长 度 区 间 分 布 的 位 置 均 发 生 了 变 
化 , 且 曲 率 变化 同步 性 较 好 的 大 长 度 区 间 的 数量 有 
所 增加 ,但 曲率 变化 同步 性 较 好 的 单一 大 长 度 区 间 
所 占 的 总 长 度 减少 了 。 这 说 明 加 入 无 迹 卡 尔 曼 滤波 
后 , 立 管 曲 率 变 化 同步 性 在 立 管 长 度 方 向 的 分 布 整 
体 性 减 小 ,离散 性 和 阶梯 性 有 所 增强 。 

总 而 言 之 ,加 入 无 迹 卡 尔 曼 滤波 后 ,整个 系统 
中 , 立 管 弯曲 变化 的 强 非 线 性 得 到 了 更 好 地 体现 。 
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时 间 /s 时 间 /s 
Et ea RP 
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寻 3(a) 不 同 管内 流 流速 下 立 管 两 端 张力 
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到 3(b) 不 同 管内 流 流速 下 立 管 张力 沿 长 度 方向 的 分 布 情况 


均匀 管内 流 流速 =0 ts 
a Mt a i 
| 均匀 管内 流 流速 =0.5 ts 


示 准 差 /kN 


立 管 张力 书 


0 200 400 600 0 200 
立 管 长 度 /m 立 管 长 度 /m 


图 3(c) 不 同 管内 流 流速 下 立 管 张力 标准 差 沿 长 度 方向 的 分 布 情况 
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图 3(d) 不 同 管内 流 流速 下 立 管 三 端 张力 谱 密 度 的 分 布 情况 
图 3 不 同 管内 流 流速 作用 下 立 管 的 张力 


上 8. 3 Dynamic response of effective tension of riser with 


XIV 


CC different velocities of uniform content flow 

有 ss 

-Cr 000 1 000 

© 均匀 管内 流 流速 =0 t/s 均匀 管内 流 流速 =0 Vs 
个 30 RiPD A 0 PD \ 
世 。- 最 小 值 | pa -最 小 值 | 
亏 09[ = 最 大 全 了 “| “名 009[3e 忆 大 信 ” | 
是 “平均 值 | | 姓 -平均 值 | 
部 400| 新 400|。 ta 
车 i | 如 CR # \ 


| 


a ) A 
OFmma ms PE | (0 mm 和 
0 200 400 600 0 200 400 600 
立 管 长 度 /m 立 管 长 度 /m 
1000 1 000 

均匀 管内 流 流速 =0.5 ts 均匀 管内 流 流速 =0.5 t/s 
合 800 
中 UKFPD 人 
名 600 “最 小 值 px 
一 最 大 值 / | 
显 4- 平 均值 了 | 
六 400 / | 
项 有 AAA 
3 200F ~ | 


0™ 350 40 60 020 40 600 
立 管 长 度 /m 立 管 长 度 /m 


到 4(a) 不 同 管内 流 流 速 下 立 管 弯 矩 沿 长 度 方向 的 分 布 情况 
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图 4(d) 不 同 管内 流 流速 下 立 管 曲率 标准 差 沿 长 度 方向 的 分 布 情况 
图 4 不 同 管内 流 流速 作用 下 立 管 的 弯 矩 和 曲率 


Fig.4 Bending moment and curvature of riser 


with different velocities of uniform content flow 

3.1.3 不 同 管内 流速 作用 下 立 管 的 转角 

图 5 为 不 同 管内 流 流速 作用 下 立 管 的 转角 变化 
情况 。 当 管内 流 为 均匀 流 时 ,对 比 观察 单一 PID 控 
制 与 UKF-PID 动力 定位 系统 作用 下 不 同 管内 流 流 
速 下 立 管 转角 RR,、R,、R, 沿 长 度 方向 的 分 布 发 现 ,加 
人 无 迹 卡 尔 曼 滤波 后 ,在 外 界 海洋 环境 载荷 下 ,滤波 
的 作用 使 得 海洋 平台 在 特定 方向 的 运动 分 量 得 到 了 
加 强 , 从 而 使 立 管 在 特定 方向 的 非 线性 弯曲 也 得 到 
了 加 强 和 体现 ,使 得 立 管 承受 绕 其 自身 y 轴 方 向 的 
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弯曲 加 强 。 这 种 现象 产生 的 原因 与 加 入 无 迹 卡 尔 曼 
滤波 后 海洋 平台 - 立 管 刚 柔 系 统 整 体 运动 的 非 线性 增 
强 有 关 , 整 体 运动 非 线性 的 增强 导致 系统 在 克服 外 界 
环境 载荷 的 变化 时 自身 的 运动 响应 变化 更 灵敏 。 
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了 5(c) 不 同 管内 流 流速 下 立 管 转角 
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四 5(f) 不 同 管内 流 流速 下 立 管 转角 R, 标准 差 沿 长 度 方向 的 分 布 情况 
色 5 不 同 管内 流 流速 作用 下 立 管 的 转角 


Fig.5 Dynamic response of rotation angle of riser with 


different velocities of uniform content flow 

对 比 观 察 单 一 PID 控制 与 UKF-PID 动力 定位 
系统 作用 下 不 同 管内 流 流速 下 立 管 转角 RR,、R, 的 标 
准 差 沿 长 度 方向 的 分 布 发 现 :加 入 滤波 后 立 管 转角 
R, 的 大 小 、 立 管 转角 R, 沿 立 管 长 度 方 向 变化 的 协 
调 性 和 同步 性 未 发 生 大 的 变化 ;加 入 无 迹 卡 尔 曼 滤 
波 后 , 立 管 转角 R, 的 变化 剧烈 程度 增加 ,但 R, 在 立 
管 长 度 方向 的 同步 性 与 协调 性 基本 没有 发 生 明 显 的 
变化 。 

对 比 观察 单一 PID 控制 与 UKF-PID 动力 定位 


系统 作用 下 不 同 管内 流 流 速 下 立 管 转角 R, 的 标准 

差 沿 长 度 方向 的 分 布 发 现 : UKF-PID 动力 定位 系统 

作用 下 的 立 管 转角 R, 变化 的 协调 性 与 同步 性 沿 立 

管 长 度 方向 没有 突变 ,但 转角 R. 变化 的 剧烈 程度 在 

长 度 方向 上 总 体 加 剧 了 。 

3.1.4 不 同 管内 流 流 速 作 用 下 海洋 平台 的 六 自由 度 
图 6 为 不 同 管内 流 作 用 下 平台 的 六 自由 度 变 
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化 。 对 比 观察 单一 PID 控制 与 UKF-PID 动力 定位 
系统 作用 下 不 同 管内 流 流速 下 的 海洋 平台 的 三 自由 
度 平 动 时 域 图 发 现 :UKF-PID 动力 定位 系统 作用 下 
的 平台 的 三 自由 度 平 动 发 生 了 较 大 的 变化 ;其 模范 
进入 稳定 状态 所 需要 的 时 间 进 一 步 缩 小 , 且 在 横 荡 
方向 达到 稳定 前 其 横 荡 的 最 大 幅度 也 大 大 减 小 ; 纵 
荡 方 向 的 波动 有 小 幅度 的 增 大 。 
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(d) 不 同 管内 流 流速 下 平台 三 自 


图 6 不 同 管内 流 流速 作用 下 平台 的 六 自由 度 


Fig. 6 Six degree of freedom of the offshore platform with different velocities of uniform content flow 
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对 比 观察 单一 PID 控制 与 UKF-PID 动力 定位 
系统 作用 下 不 同 管内 流 流速 下 的 平台 的 三 自由 度 转 
动 时 域 图 像 发 现 ，UKF-PID 动力 定位 系统 作用 下 的 
平台 三 自由 度 转动 的 时 域 图 的 主要 变化 集中 体现 在 
赠 摇 角 的 变化 上 ,平台 的 秋 摇 角 在 时 域 上 经 过 最 初 
的 时 间 步 上 的 短暂 停滞 后 , 在 后 续 的 时 间 步 上 有 了 
小 幅度 的 持续 波动 。 产 生 这 种 现象 的 原因 也 与 加 入 
无 迹 卡 尔 曼 滤波 后 系统 整体 的 非 线性 增强 有 关 。 
对 比 观察 单一 PID 控制 与 UKF-PID 动力 定位 
系统 作用 下 不 同 管内 流 流速 下 的 海洋 平台 的 三 自由 
度 平 动 谱 密度 图 像 可 以 发 现 :平台 平 动 谱 密度 的 变 
化 主要 集中 在 纵 荡 方 向 ,但 在 UKF-PID 动力 定位 系 
统 作 用 下 ,海洋 平台 的 纵 荡 谱 密度 只 有 一 个 极 大 的 
高 峰值 和 一 个 很 小 的 低 峰值 ;也 就 是 说 ,在 UKF-PID 
动 为 定位 系统 作用 下 ,虽然 纵 荡 谱 密度 的 峰值 有 所 
出 加 ,但 平台 纵 荡 运 动 可 能 出 现 的 超 谐 共振 消失 了 ， 
这 对 于 提高 整个 海洋 平台 - 立 管 刚 柔 系统 的 安全 性 
是 有 益 的 ,同时 也 降低 了 疲劳 损伤 发 生 的 可 能 | 
3 段 塞 流 时 的 计算 结果 
在 本 小 节 的 仿真 中 ,具体 立 管内 的 段 塞 流 参数 


着 役 塞 流 流速 为 0.5 m/s, 段 塞 点 之 间 的 流体 密度 
海 8135 tm , 立 管 上 出 现 段 塞 流 的 长 度 为 320 m, 段 
号 站 处 的 总 密度 为 0. 5 wm ,每 两 个 相 邻 的 段 塞 点 
之 加 的 距离 为 10 m, 段 塞 点 的 个 数 为 52 个 ; 段 塞 流 
的 赵 始 计算 参考 点 设 为 长 度 为 11 m 的 滑 移 接 头 端 ， 
渭 鼻 接 头 端 是 没有 段 寒流 的 。 

3.21 段 塞 流 作用 下 立 管 的 张力 

三 图 7 为 段 寒流 作用 下 的 立 管 张力 变化 情况 。 对 
比 岗 察 单一 PID 控制 与 UKF-PID 动力 定位 系统 作 
用 时 段 塞 流 作用 下 , 立 管 的 张力 大 小 ` 张 力 标准 差 及 
张力 谱 密 度 的 分 布 情况 可 以 发 现 :两 种 控制 模式 下 
均匀 管内 流 作用 下 的 计算 结果 相差 不 大 ; 与 单一 
PID 动力 定位 系统 作用 时 段 塞 流 作用 下 的 对 应 计算 
结果 相 比 ,UKF-PID 控制 模式 下 的 立 管 张力 大 小 和 


时 间 /s 


7(a) 段 塞 流 时 立 管 两 端 张力 
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到 7(b) 段 塞 流 时 立 管 张力 沿 长 度 方 向 的 分 布 情况 
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图 7 段 塞 流 作用 下 立 管 的 张力 


Fig.7 Dynamic response of effective tension of riser 


with different velocities of slug flow 

3.2.2 段 塞 流 作 用 下 立 管 的 谊 矩 和 曲率 
图 8 为 段 塞 流 作用 下 立 管 的 弯 矩 和 曲率 分 布 情 
况 。 由 图 8 可 知 ,在 UKF-PID 动力 定位 系统 作用 下 ， 
立 管 弯 矩 和 曲率 沿 立 管 长 度 方向 的 分 布 情况 与 
UKF-PID 动力 定位 系统 作用 时 均匀 管内 流 流速 下 立 
管 的 弯 矩 和 曲率 沿 立 管 长 度 方向 的 分 布 情况 相 比 变 
化 较 大 ,这 种 变化 主要 集中 于 弯 和 矩 和 曲率 的 大 小 变 
化 上 ,此 时 同一 位 置 的 立 管 弯 矩 最 大 值 . 曲 率 最 大 值 
普遍 减 小 ,但 曲线 的 总 体形 态 与 均匀 管内 流 作 用 时 
相差 不 大 。 但 同 单一 PID 动力 定位 系统 作用 时 段 塞 
流 作 用 下 立 管 弯 矩 和 曲率 治 立 管 长 度 方向 的 分 布 情 
况 相 比 , UKF-PID 动力 定位 系统 作用 时 段 塞 流 作用 
下 同一 位 置 的 立 管 弯 矩 最 大 值 、 曲 率 最 大 值 普遍 增 
大 , 且 弯 和 矩 和 曲率 在 立 管 长 度 上 的 分 布 曲线 的 形态 
也 发 生 了 一 定 程度 的 变化 ,虽然 这 种 变化 并 不 大 ,但 
在 整个 长 度 范围 内 , 立 管 承受 最 大 或 最 小 这 曲 作用 


的 部 位 还 是 发 生 位 置 前 移 或 后 移 , 在 实际 工程 中 这 


种 距离 偏 移 还 是 不 能 忽视 的 。 
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(d) 不 同 管内 流 流速 下 立 管 曲率 标准 差 沿 长 度 方向 的 分 布 情况 
8 段 塞 流 作 用 下 立 管 的 弯 矩 和 曲率 
Fig.8 Bending moment and curvature of riser 


with different velocities of slug flow 
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3.2.3 段 塞 流 作 用 下 立 管 的 转角 

图 9 为 段 塞 流 作 用 下 立 管 的 转角 变化 情况 。 由 
图 9 可 知 :与 单一 PID 动力 定位 系统 作用 时 段 塞 流 
作用 下 立 管 转角 R,、R,、R. 相 比 ，UKF-PID 动力 定位 
系统 作用 下 的 转角 R, 大 小 沿 立 管 长 度 方向 的 分 布 
情况 基本 不 变 ;转角 RR, 的 最 大 值 变 大 ,最 小 值 变 小 ， 
且 转 角 R, 沿 立 管 长 度 方向 的 曲线 形态 也 发 生 了 明 
显 的 变化 ;在 UKF-PID 动力 定位 系统 作用 下 , 立 管 
同一 位 置 ,转角 R,、R, 数值 上 相差 增 大 , 且 转 角 尺 ， 
明显 大 于 转角 RR,, 这 说 明 在 此 时 立 管 发 生 的 弯曲 主 
要 以 绕 y 轴 的 弯曲 为 主 。 

与 单一 PID 动力 定位 系统 作用 时 段 塞 流 作 用 下 
时 的 情况 相 比 , 立 管 转角 R,、R, 变化 在 立 管 长 度 方 
向 不 存在 同步 性 和 协调 性 ; 且 由 于 此 时 同一 位 置 转 
角 R,.R, 标准 差 数值 上 也 相差 较 大 ,因此 此 时 立 管 
同一 位 置 的 转角 R, 变化 和 转角 R, 变化 也 不 存在 协 
调 性 与 同步 性 。 这 种 立 管 同一 位 置 转角 R,、R, 变化 
的 非 同 步 性 与 非 协调 性 某 种 程度 上 来 说 ,会 Sm. 
些 部 位 的 磨损 ,增加 这 些 部 位 的 寿命 ,从 而 减少 立 
的 疲劳 损伤 。 而 UKF-PID 和 单一 PID 动力 定位 系 
统 作用 时 段 塞 流 作 用 下 立 管 转角 R. 标准 差 沿 立 管 
长 度 方 向 的 分 布 差别 不 大 ,这 说 明 , 当 管内 流 形式 为 
段 塞 流 时 ,两 种 控制 方式 下 立 管 转角 RR, 变化 的 同步 
性 相差 不 大 ; 沿 立 管 长 度 方向 立 管 绕 自身 轴 向 的 扭 
转 都 能 保持 较 好 的 同步 性 和 协调 性 ,这 种 立 管 整 体 
长 度 方向 扭转 变化 同步 性 对 于 保证 立 管 稳定 的 在 位 
运行 是 有 积极 作用 的 。 

0.8 
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外 RR, 沿 长 度 方向 的 分 布 情况 
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9(a) 段 塞 流 时 立 管 转 


图 8 中 ,UKF-PID 控制 模式 下 沿 立 管 长 度 方向 
弯 矩 标准 差 数值 上 明显 增 大 ,这 说 明 对 于 立 管 的 具 
体 某 一 位 置 而 言 ,其 弯曲 变化 在 这 种 状态 下 更 加 急 
剧 ; 但 对 于 整个 立 管 长 度 方向 上 的 弯曲 变化 来 说 ,由 
于 除 项 端 和 尾 端 外 , 立 管 任 一 位 置 的 弯曲 变化 的 急 
剧 程度 都 增加 了 ,因而 其 同步 性 和 协调 性 沿 立 管 整 
体 长 度 方向 对 应 的 同步 性 和 协调 性 相差 并 不 大 。 
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9(b) 段 塞 流 时 立 管 转 f 
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立 管 长 度 /m 
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了 9(c) 段 塞 流 时 立 管 转角 R. 沿 长 度 方向 的 分 布 情况 
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3 段 宏 流 作用 下 的 海洋 平台 六 自由 度 
CN 图 10 为 段 赛 流 作用 下 海洋 平台 的 六 自由 度 
人 人 这 。 对 比 观察 图 10 中 海洋 平台 的 六 自由 度 运 
动 时 域 响应 计算 结果 可 以 发 现 :在 时 域 上 ,平台 的 模 
RE 
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提前 ,平台 纵 荡 幅度 也 有 所 增加 ,而 其 身 播 角 在 
时 右上 也 出 现 了 极 小 幅度 的 有 序 持 续 性 波动 。 
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10(a) 段 塞 流 时 平台 三 自由 度 平 动 时 域 曲线 
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图 10(b) 段 塞 流 时 平台 三 自由 度 平 动 谱 密 度 
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到 10(d) 段 塞 流 时 平台 三 自由 度 转动 谱 密 度 

图 10 段 塞 流 作用 下 的 海洋 平台 六 自由 度 
Fig. 10 Six degrees of freedom of the offshore 


platform under slug flow 

观察 图 10 中 海洋 平台 的 三 自由 度 平 动 谱 密 度 
图 像 可 以 发 现 ,在 UKF-PID 动力 定位 系统 作用 下 ， 
虽然 纵 荡 谱 密度 的 峰值 有 所 增加 ,但 平台 纵 荡 运 动 
的 超 谐 共振 消失 了 (这 一 结论 在 管内 流 形式 为 均匀 
流 时 也 是 成 立 的 ) ,这 对 于 提高 整个 海洋 平台 - 立 管 
刚 柔 系 统 的 安全 性 是 有 益 的 ,同时 也 降低 了 疲劳 损 
伤 发 生 的 可 能 性 。 也 就 是 说 ,在 这 种 状态 下 ,海洋 平 
台 通 过 三 自由 度 平 动 运动 释放 外 界 载 荷 能 量 的 过 程 
中 ,通过 纵 荡 释放 掉 的 能 量 所 占 的 比例 增加 ,通过 垂 
荡 运 动 释放 掉 的 能 量 所 占 的 比例 有 所 减少 。 

对 比 观 察 单一 PID 控制 和 UKF-PID 动力 定位 
系统 作用 时 分 别 在 段 塞 流 作 用 下 与 均匀 管内 流 流速 
作用 下 海洋 平台 三 自由 度 谱 密 度 图 像 可 发 现 ,在 段 
塞 流 作 用 下 ,两 种 控制 模式 下 海洋 平台 的 纵 荡 谱 密 
度 峰 值 均 较 均匀 管内 流 时 有 所 增加 。 这 是 因为 管内 
流 形式 的 改变 造成 了 整个 海洋 平台 - 立 管 系统 的 总 
体 刚 度 发 生变 化 ,进而 使 在 运动 非 线 性 增强 的 前 提 
下 海洋 平台 的 纵 荡 方 向 的 约束 减弱 ,这 也 证 明 管 内 
流 的 形式 对 整个 系统 的 响应 存在 影响 。 


4 结 论 


刚 柔 多 体系 统 的 相互 作用 是 一 个 非常 复杂 的 、 
不 可 预测 的 、 非 线性 的 作用 。 随 着 立 管 中 管 内 流 流 
速 和 管内 流 形 式 的 变化 , 立 管 挠 性 的 变化 使 立 管 与 


海洋 平台 之 间 的 相互 作用 发 生变 化 ,从 而 对 海洋 平 
台 的 动力 定位 过 程 造成 影响 ,最 终 影响 整个 动力 定 
位 海洋 平台 - 立 管 系统 的 总 体 响应 。 

动力 定位 系统 的 形式 .管内 流 的 形式 都 会 对 立 
管 弯 曲 的 传递 情况 和 沿 立 管 长 度 方向 变化 的 同步 
性 .协调 性 及 波动 的 急剧 性 产生 重要 影响 ,其 中 管内 
流 的 形式 起 着 主要 作用 ,动力 定位 系统 的 作用 也 不 
可 忽视 。 

引起 平台 三 自由 度 平 动 谱 密度 发 生变 化 的 主要 
原因 是 定位 控制 中 加 入 了 无 迹 卡尔 曼 滤波 ,这 激发 
了 海洋 平台 三 自由 度 平 动 运动 非 线性 的 增强 ,并 改 
变 了 能 量 在 海洋 平台 三 种 平 动 运动 的 分 布 情况 。 

加 入 无 迹 卡 尔 曼 滤波 后 ,海洋 平台 - 立 管 刚 柔 系 
统 整 体 上 在 克服 外 界 环境 载荷 的 变化 时 的 运动 响应 
变化 更 灵敏 。 
[一 由 于 动力 定位 系统 在 定位 过 程 中 一 般 只 考虑 平 
全 浊 动 的 低频 部 分 ,因此 平台 的 高 频 运动 信号 以 及 
答 种 干扰 信号 需要 通过 无 迹 卡尔 曼 滤 波 消除 ; 除 此 
区 天 , 由 于 无 迹 卡尔 曼 滤波 自身 的 功能 作用 ,也 可 以 
使 得 整个 海洋 平台 - 立 管 系统 的 非 线性 得 到 更 好 地 
依 现 ,避免 某 种 程度 上 的 强 非 线性 结果 线性 化 。 
> 
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